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Interchain coupling effects on photoinduced phase transitions between neutral and 
ionic phases in an extended Hubbard model with alternating potentials and an 

electron-lattice coupling 
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Dynamics of ionic-to-neutral and neutral-to-ionic phase transitions induced by intrachain charge- 
transfer photoexcitations are studied in a quasi-one-dimensional extended Hubbard model with 
alternating potentials and an electron-lattice coupling for mixed-stack charge-transfer complexes. 
For interchain couplings, we use electron-electron interactions previously estimated for TTF-CA 
(TTF=tetrathiafulvalene, CA=chloranil) . Photoexcitation is introduced by a pulse of oscillating 
electric field. The time-dependent Hartree-Fock approximation is used for the electronic part, and 
the classical approximation for the lattice part. In the ionic-to- neutral transition, the transferred 
charge density is a strongly nonlinear function of the photoexcitation density, which is characterized 
by the presence of a threshold. With substantial interchain couplings comparable to those in TTF- 
CA, the interchain correlation is strong during the transition. Neutral domains in nearby chains 
simultaneously grow even if their nucleation is delayed by reducing the amplitude of the electric field. 
With weaker interchain couplings, the growing processes are in phase only when the amplitude of 
the electric field is large. Thus, the experimentally observed, coherent motion of a macroscopic 
neutral-ionic domain boundary is allowed to emerge by such substantial interchain couplings. In the 
neutral-to-ionic transition, by contrast, the transferred charge density is almost a linear function 
of the photoexcitation density. Interchain electron-electron interactions make the function slightly 
nonlinear, but the uncooperative situation is almost unchanged and consistent with the experimental 
findings. 

PACS numbers: 71.10.Hf, 71.35.-y, 63.20.Kr, 78.47.-|-p 
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I. INTRODUCTION 

Photoinduced phase transitions have attracted much 
attention i Initially local structural deformations trigger 
a macroscopic change in dielectric, optical, and/or mag- 
netic properties. The transition dynamics are recently 
studied in different molecular materials, revealing the ori- 
gin of cooperativity responsible for the proliferation of 
induced electronic states.^ These researches would lead 
to developing techniques for dynamically controlling the 
electronic and structural properties and their coherence 
in such materials. 

The mixed-stack organic charge-transfer complex, 
TTF-CA, is one of the most intensively studied mate- 
rials among them. The donor TTF and acceptor CA 
molecules are alternately stacked along the most con- 
ducting axis. At low temperature or under high pressure 
with contraction, they are ionic due to the long-range 
Coulomb interactioni^ Otherwise, they are neutral due 
to the difference between the ionization potential of the 
donor molecule and the electron affinity of the acceptor 
molecule. In the ionic phase at ambient pressure, these 
molecules are dimerized. 

Both ionic-to-ncutral and ncutral-to-ionic transitions 
are induced by photoirradiation near the energy of the lo- 
calized intramolecular excited state of TTF moleculesi^ 
In the latter transition, the photoinduced phase is re- 
cently shown by X-ray diffraction to be accompanied with 
three-dimensionally ordered ferroelectric polarizationsi^ 



The ionic-to-neutral transition induced by intrachain 
charge-transfer photoexcitations has a threshold in the 
photoexcitation density, below which a macroscopic neu- 
tral domain cannot be generated.^ 

Ultrafast optical switching from the ionic to the neu- 
tral states is observed after a resonant excitation of 
the charge-transfer band polarized along the stacking 
axisi Recently, experimental results have been summa- 
rized covering both the ionic-to-neutral and the neutral- 
to-ionic transitions with different excitation energies, 
excitation densities, and temperatures^ In particular, 
the dynamics of the neutral-to-ionic transition induced 
by charge-transfer photoexcitations polarized along the 
stacking axis is clearly different from that of the ionic- 
to-neutral one. Although ionic domains arc initially pro- 
duced by lights, they quickly decay even if the excitation 
density is high. The initial conversion fraction is a linear 
function of the excitation density. 

The photoinduced phase transitions in TTF-CA have 
been studied also theoretically,^'^" and the origin of the 
different dynamics in the two transitions is discussed . ^^i — 
However, these studies are limited to one-dimensional 
electron models. Although the short-time behavior may 
not be so sensitive to interchain couplings, they would 
eventually alter the long-time dynamics. For instance, 
they are crucial to the emergence of coherence, which can 
be manifested by introducing a double pulse .--'-— ' The ef- 
fect of interchain elastic couplings was considerediSi How- 
ever, interchain electron-electron interactions are much 



stronger and would be much more influential. The ab ini- 
tio quantum chemical calculation for the intramolecular 
charge distribution in TTF and CA moleculesi^ suggests 
that the interchain electrostatic energies between neigh- 
boring molecules are smaller than but comparable to the 
intrachain ones. 

In this paper, we add interchain electron-electron in- 
teractions to the previous model for TTF-CA and study 
their effects on the phase transition dynamics induced by 
intrachain charge-transfer photoexcitations. The qualita- 
tive difference found in the purely one-dimensional model 
is basically unchanged: the ionic-to-neutral transition 
proceeds in a cooperative manner, whereas the neutral- 
to-ionic one in an uncooperative manner. However, the 
dynamics of the ionic-to-neutral transition strongly de- 
pend on the strength of interchain couplings. Indeed, 
substantially strong interchain electron-electron interac- 
tions are important to maintain the coherent motion of 
the neutral-ionic domain boundary. They allow neutral 
domains in nearby chains to grow simultaneously even if 
their nucleation is delayed by reducing the amplitude of 
the electric field. 



II. QUASI-ONE-DIMENSIONAL MODEL FOR 
MIXED-STACK CT COMPLEXES 



For electrons in the highest occupied molecular or- 
bital (HOMO) at donor sites and the lowest unoccupied 
molecular orbital (LUMO) at acceptor sites, we use a 
quasi-one-dimensional extended Hubbard model with al- 
ternating potentials and an electron-lattice coupling at 
half filling^ 
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and with the interchain component, 



+Vp2Snij6ni+ 



ij+i > 



(2) 



(3) 



where, cj ■ ^ {(^id.cr) is the creation (annihilation) opera- 
tor of an electron with spin a at site I of chain j , nij_cF = 
C;j,^Qj,o., nij = n/j^i -I- n/j^|, uij is the dimension- 
less lattice displacement of the molecule from its equidis- 
tant position along the chain, and yij = u/-f-ij — uij. 
The distance between the Ith and (^+l)th molecules is 
then given by rij = ro(l -I- ui^ij — ui,j) along the jth 
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FIG. 1: Schematic electronic and lattice structures in (a) neu- 
tral and (b) ionic phases. 



chain, where rg is the averaged distance between the 
neighboring molecules. The intersite Coulomb interac- 
tions are not between the electron densities, nij, but 
between the excess electron densities, 5nij, which we de- 
fine as 5ni,.j = ni_j — 2 for odd I (at a donor site) and 
5nij = nij for even I (at an acceptor site). The to- 
tal charge of the donor molecule at site I is then given 
by —eSnij — +e{2 — n-ij), while that of the acceptor 
molecule by —eSnij = —enij. 

The nearest-neighbor interaction strength along the 



chain Vij depends on the bond length yij, Vi^ 



V 



l32yij, where V is for the regular lattice. The param- 
eter ig denotes the nearest-neighbor transfer integral, d 
the level difference between the neighboring orbitals in 
the neutral limit, and U the on-site repulsion strength. 
The elastic energy is expanded up to the fourth order 
with the parameters fci and ^2- The quantity mij de- 
notes the molecular mass. The interchain interactions 
are characterized by the repulsion strength Up for neigh- 
boring donor molecules and for neighboring acceptor 
molecules, I^pi and V^2 for neighboring donor and ac- 
ceptor molecules. Because of the molecular tilting, Vpi 
and V^2 are not equal. The number of sites in a periodic 
chain is denoted by L, and that of chains by M. 

In the neutral phase, the orbital of the donor molecule 
is nearly doubly occupied, while that of the acceptor 
molecule is nearly empty [Fig. ^a)] . The total charge 
of any molecule is close to zero. In the ionic phase, 
both orbitals are nearly singly occupied [Fig.QJb)]. The 
total charges of the donor and acceptor molecules are 
close to ±e. The ionicity is defined as p = 1 + 

l/(LM)EtiE-ii(-l)'K.)- 



III. EFFECTS OF INTERCHAIN COUPLINGS 

In this section only, the effects of the interchain cou- 
plings are discussed for io=/32=0. The edge chains j=l 
and j=M are ignored unless otherwise stated. The total 
energy per donor-acceptor unit is Ejsi = —d + U in the 
neutral phase [Fig.UJa)], and Ei = ~2V + 2Up - 2Vpx - 
2Vp2 in the ionic phase [Fig.^b)]. The energy difference 
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FIG. 2: Neutral domain in ionic background. 



is then given by 

En-Ei = U + 2V- rfren , 

with 

dron = d+2Up- 2ypi - 2yp2 



(4) 



(5) 



The ab initio quantum chemical calculation has es- 
timated the electrostatic energies between neighboring 
molecules as F=1.26 eV, f/p=0.94 eV, and Vpi = 1.08 eV 
for TTF-CA.J^ The inequahty 1/pi > Up, i.e., the domi- 
nance of the interchain attraction over the interchain re- 
pulsion in the ionic phase, is responsibleii for the phase 
diagram containing the paraelectric ionic phasei^ and for 
the discontinuous change of the lattice parameter b at 
the neutral-ionic transitionriS The inequality is due to 
the molecular tilting, so that T/p2 is much smaller than 
Ipi. Because the above estimations are based on the in- 
tramolecular charge distribution for isolated molecules, 
the estimated values should be regarded as upper bounds 
for the interaction strengths. In numerical calculations, 
we will fix t/p/Vpi=0.9 and Vp2=0 and vary Up and Vpi. 
Here, we do not fix them for general discussions. 

To create a neutral domain of Ld donor-acceptor units 
in the ionic background (Fig. [SJ , it costs 



V + L42Vpi + 2Vp2 ~ 2Up) 



(6) 



at the phase boundary, E^=Ei. The domain wall energy 
is then given by E'dw = V/2^ as in the one-dimensional 
caseiSi The confinement of a metastable domain is given 
by the id-dependent term, in which the large contribu- 
tion from the Up term nearly cancels out the largest one 
from the V^i term due to the relation, V^i ^ Up :s> Vp2. 
The fact that they are comparable is important to facil- 
itate the growth of metastable domains. 

The optical gap for the excitation along the stacking 
axis is given by 



ET = drcn -U-V-2Up 
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in the neutral phase [Fig. ISJa)], and by 
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in the ionic phase [Fig.jSJb)]. Here also, the optical gaps 
are not drastically modified by the interchain coupHngs 
owing to the relation, T/pi ^ Up ^ Fp2. For drcn, we 
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FIG. 3: Intrachain charge-transfer excitations in (a) neutral 
and (b) ionic phases. 



will employ the value that was previously used for d in 
the one-dimensional case. Otherwise, we do not need 
to modify the model parameters estimated previously. 
Now, we consider the edge chains j=l and j=M. The 
above equation for the optical gap is not modified in the 
neutral phase. In the ionic phase, the optical gap at the 
edge chains is given by 



E] 



lOpt 



+ U + 3V , 



(9) 



which is slightly smaller than the bulk contribution. 



IV. METHOD OF CALCULATING DYNAMICS 

The method is the same as that used previouslyfii*i2ii^ 
including the time-dependent Schrodinger equation for 
the electronic part within the unrestricted HF approx- 
imation, the classical equation of motion for the lat- 
tice part according to the Hellmann-Feynman theorem, 
the addition of random numbers to the initial yij and 
ill J values obeying the Boltzmann distribution at a fic- 
titious temperature T, and the incorporation of a time- 
dependent electric field. 



E{t) = ii^cxtsinwcxti , 



(10) 



with amplitude ii^oxt and frequency Wcxt for < i < 
2nNeict/^e-Kt with integer A^oxt, into the Peierls phase of 
the transfer integral. E{t) is zero otherwise. 

Because the infinitesimal deviations from the static 
self-consistent HF solution are equivalent to individ- 
ual/collective excitations coupled with phonons in the 
random phase approximation (RPA),'= the present time- 
dependent HF calculations treating finite deviations nat- 
urally go beyond the RPA including finite lifetimes of 
excitations.^^ The linear absorption spectra have peaks 



at the energies of excitons in the RPA, which are broad- 
ened owing to interactions among excitons and other ex- 
citations. The RPA itself is beyond the Tamm-DancofF 
approximation, which is nothing but the single-excitation 
configuration-interaction method. In the previous sec- 
tion, the effects of interchain couphngs are discussed in 
the strong-couphng hmit. They are described in the 
Tamm-Dancoff approximation also. Indeed, the present 
HF method produces the linear absorption spectra that 
are consistent with these strong-coupling estimations: 
it is not really a weak-coupling approach but a small- 
amplitude-fluctuation approach. It goes beyond the 
rigid-band picture by containing dynamical fluctuations 
but cannot qualitatively describe quantum spin fluctua- 
tions after artificial breaking of the spin rotational sym- 
metry. Because the relevant energy scale is large in 
the photoinduced processes, the poor approximation for 
spin fluctuations would not qualitatively alter the charge- 
lattice coupled dynamics of present concern. 

In the one-dimensional case, the model parameters 
are chosen to reproduce the ab initio estimation of the 
transfer integral and the measured values of the ionic- 
ity, the dimerization, and the absorption spectral We 
employ basically the same parameters as beforeiSiiiiiSii^ 
for M=10 chains of L=100 sites: io=0.17 eV, C/=1.528 
eV, F=0.604 eV (when the ionic phase is photoex- 
cited) or F=0.600 eV (when the neutral phase is pho- 
toexcited), d,c^=2.716 eV /32=8.54 eV, fci=4.86 eV, 
^2=3400 eV, and the bare phonon energy ujopt = 
(l/ro)(2A:i/mr)i/2=0.0192 eV-M. Here, the reduced mass 
rrir is defined as mr = miymx/ {rriD + toa) with tod for 
the donor molecule and ttia for the acceptor molecule. 
Note that the phonon energy used here is a couple of 
times larger than that in TTF-CA. The main results are 
not altered by the choice of its value. 

If we adopt the ratios U-p/V and V-pi/V evaluated by 
the ah initio quantum chemical calculationi^ and the 
value y=0.604 eV above for the ionic phase, the inter- 
chain couplings are roughly estimated as t/p=0.45 eV and 
Vpi=0.50 eV. Hereafter these parameters are written in 
units of eV. The ionicity in the static self-consistent HF 
solution is 0.20 (0.18) for Fpi^O (Vpi=0.50) in the neu- 
tral phase and 0.95 (0.97) for Vpi=0 (ypi=0.50) in the 
ionic phase, which are closer to the strong-coupling lim- 
its than the experimentally estimated values of about 0.3 
and 0.7. 



V. RESULTS 
A. lonic-to-neutral transition 

The real-time dynamics shown later depend on the fre- 
quency of the electric field. So, the linear absorption 
spectra in the ionic phase are calculated for reference. 
Without interchain couplings, the position of the absorp- 
tion peak is located at Wext ~ 32ci;opt ~0.6I eV.^-^ With 
interchain couplings t/p=0.45 eV, T/pi=0.5 eV and Vp2=0, 



it would be shifted upward by (2Vpi -I- 2Vp2 — 2Up)=Q.\ 
eV'^ Swopt according to Eq. JSJ in the strong-coupling 
limit. The highest peak is actually located at Wcxt ^ 
38ct;opt, corresponding to the charge-transfer excitons in 
the bulk. The second highest peak is at Wext ~ 32ct;opt 
owing to the charge-transfer excitons at the edges. The 
latter position is the same with the one-dimensional case, 
as Eq. (PJ suggests. These assignments are confirmed by 
the fact that resonant photoexcitations produce neutral 
domains in the respective positions. 

When the ionic phase is photoexcited, the ionicity 
(O/j, defined as pij == 1 + (-I)'(-(n/_ij) -f 2(n;j) - 
('T-i+ij})/4, changes abruptly from a large value to a 
small value during the transitionJi Of course the tim- 
ing of this change depends on the place. Thus, it is easy 
to trace the development of neutral domains. After a 
neutral domain is created in a chain, it may grow to con- 
vert the whole chain into the neutral phase or may be 
suppressed by nearby ionic chains to return to the ionic 
phase. For each chain, the final state is either neutral 
or ionic. However, it is possible for the final state to 
be a mixture of neutral and ionic chains if the supplied 
energy from the oscillating electric field is insufficient. 
Below we plot the number of chains converted from the 
ionic phase to the neutral phase among ten chains. Once 
a chain is converted, further application of the electric 
field may increase the ionicity even when the transition 
is not completed in other chains. This is not really a 
photoinduced effect, but rather a thermal effect in that 
the excess energy randomizes the electron distribution. — 
It indeed happens for isolated or weakly coupled chains. 
Therefore, we have carefully looked into the evolution of 
ionicity at each chain. 

The number of photoconverted chains is plotted in 
Fig. ^ as a function of the number of absorbed pho- 
tons per chain, which is obtained by dividing the incre- 
ment of the total energy per lOO-site chain by the fre- 
quency of the electric field. The frequency of the electric 
field is so varied with interchain couplings, Vpi=0.3, 0.4, 
and 0.5 with C/p/V^i=0.9, that the energy of photons 
relative to that of the linear absorption peak is almost 
unchanged. To convert only one chain among ten, the 
strongly coupled chains need the higher density of pho- 
tons. This is reasonable because the strongly coupled 
chains need a higher energy to create a neutral domain 
in the ionic background, according to Eq. ©. To con- 
vert the whole system, however, the strongly coupled 
chains need a much lower density of photons than the 
weakly coupled chains. Once a neutral domain is nucle- 
ated, other neutral domains are almost simultaneously 
created in the neighboring chains to grow spontaneously. 
Their growth cannot be stopped even if the electric field 
is switched off immediately after the appearance of the 
first domain. As a consequence, the final state in the 
strong-coupling case is either a globally neutral one or a 
globally ionic one. That is why the corresponding plot is 
a step function. 

In order to see how the mixture of neutral and ionic 
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FIG. 4: (Color online) Number of converted chains, as a func- 
tion of number of absorbed photons per chain, for different 
Vpi and Ci;cxt/<^opt- The electric field with e_Eoxt'"o/'^opt=6 is 
applied to the ionic phase at T/to—10~^. 
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FIG. 5: (Color online) Number of converted chains, as a func- 
tion of number of absorbed photons per chain, for isolated 
chains, Vpi=0, with a;ext/<^opt=31. The electric fields with 
eEcxtro/i^opt=4:, 5, and 6 are applied to the ionic phase at 
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chains are produced for weak interchain interactions, we 
uncouple the chains by setting Vpi=0 and use the same 
random numbers for the initial yi^ and iiij values as in 
the cases of finite interchain interactions. The number 
of photoconverted chains is plotted in Fig. |S1 Because 
all the chains have the same condition except for the 
weakly randomized, initial lattice variables, the mixture 
of neutral and ionic chains is due to the distribution of 
the threshold photoexcitation density. Photoexcitation 
densities in some chains are above the threshold, while 
those in other chains are below the threshold. The dis- 
tribution is substantially broad for small amplitude of 
the electric field. Within a chain, the transition proceeds 
cooperatively once a sufficiently long neutral domain is 
nucleated. The nucleation processes are randomly trig- 



gered by lattice fluctuations, so that they occur at differ- 
ent times. Different chains have absorbed different den- 
sity of photons before the nucleation process. The distri- 
bution of the threshold photoexcitation density becomes 
narrower as the amplitude of the electric field increases. 
In Fig. [SI some plots are not monotonic, where an ini- 
tially created neutral domain happens to be completely 
annihilated by continuing photoirradiation and another 
neutral domain is created later in the same chain but at a 
different place, which eventually grows and finally covers 
the whole chain lii 

For weak interchain electron-electron interactions, the 
relation between the number of photoconverted chains 
and the number of absorbed photons is similar to that for 
isolated chains (not shown). Indeed, for small amplitude 
of the electric field, the range of the photoexcitation den- 
sity for producing a mixture of neutral and ionic chains 
is wide. As the amplitude of the electric field increases, 
this range becomes narrower. The interchain correlation 
is very weak during the transition as long as the elec- 
tric field is not so strong. However, the way in which 
the transition proceeds is different from that for isolated 
chains. When the frequency of the electric field is closer 
to the energy of charge-transfer excitons in the bulk than 
that at the edges, a first neutral domain is created around 
the central chains by the low density of photons. If the 
electric field is switched off immediately after the first do- 
main appears, only the chain containing this domain is 
converted into a neutral one, and the residual chains re- 
main ionic. By continuing the application of the electric 
field, next domains are created in the neighboring chains. 
The positions of their appearances are random along the 
chains for weak interchain electron-electron interactions. 

Comparing the results for Vpi=0.3 and those for 
Vpi—QAj^^ we find that the interchain electron-electron 
interactions increase both the probability of suppression 
after the appearance of a first domain by neighboring 
ionic chains and the probability of nucleation of a second 
domain near the first one in the neighboring chains. The 
interchain correlation is accordingly increased during the 
transition dynamics. Nevertheless, the interchain corre- 
lation is generally very weak unless the amplitude of the 
electric field is large. For very large amplitude of the elec- 
tric field, metastable domains appear quickly and almost 
simultaneously after the photoirradiation starts even if 
the interchain electron-electron interactions are weak. 
Such coherence that is forced by the intense electric field 
survives shortly and decays rather quickly. This situation 
may be related to the coherent oscillation observed in the 
photoexcited quarter-fiUed-band charge-ordered organic 
(EDO-TTF)2PF6 saltiSa 

For strong interchain electron-electron interactions 
comparable to those in TTF-CA, the number of photo- 
converted chains is always a step function of the number 
of absorbed photons (Fig. 0) . The threshold increment of 
the total energy slightly increases with decreasing ampli- 
tude of the electric field. For small amplitude of the field, 
it takes a very long time to absorb the sufficient amount 
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FIG. 6: (Color online) Number of converted chains, as a func- 
tion of number of absorbed photons per chain, for strongly 
coupled chains, Vpi=0.5, with a;ext/<^opt=36. The electric The electric field with e-Bext?'o/'^opt=3 is applied to the neu- 
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FIG. 7: (Color online) Final ionicity, as a function of number 
of absorbed photons per chain, for different Vpi and Wext/wopt- 



fields with e-Eext'"o/'^opt=3, 4, and 6 are applied to the ionic 
phase at r/io=10"^. 



tral phase at T/to- 
chains.^'^ 



10" 



The straight line is for isolated 



of energy to nucleate a neutral domain. Even if the ap- 
pearance of a first neutral domain is delayed by such re- 
duction of the field amplitude, it is promptly followed by 
the appearances of other neutral domains near the first 
one in the neighboring chains. It is like an avalanche 
phenomenon. The growth of these domains is not sup- 
pressed by the surroundings any more. The growth of 
neutral domains coherently proceeds, which is consistent 
with the experimentally observed, coherent motion of the 
macroscopic neutral-ionic domain boundary on a large 
time scaled 



B. Neutral-to-ionic transition 

The linear absorption spectra in the neutral phase are 
first calculated before discussing the photoinduced dy- 
namics from the neutral phase. Without interchain cou- 
plings, the position of the absorption peak is located at 
(^cxt ~ 30wopt ~0.58 eVJ^^ With interchain couplings 
[/p=0.45 eV, ypi=0.5 eV and Vp2=0' it would be shifted 
upward again by (2Vpi -I- 2Vp2 — 2t/p)=0.1 eV~ 5wopt 
according to Eq. iQ in the strong-coupling limit. The 



peak is actually located at cjext ~ 34woptj corresponding 
to the charge-transfer excitons both in the bulk and at 
the edges. 

Then, the neutral phase is photoexcited. In contrast 
to the ionic-to-neutral transition, the ionicity changes 
gradually from a small value to a large value during 
the neutral-to-ionic transition. After the electric field 
is switched off, photoinduced ionic domains spatially 
fluctuatei^ and the averaged ionicity slightly relaxes, but 
no large variation is observed in charge and lattice dy- 
namics even if the transition is not completed. In other 
words, the photoinduced ionic domains do not grow with- 
out energy supply from the oscillating electric field. It 



implies that the neutral-to-ionic transition proceeds in an 
uncooperative manner. We plot the ionicity in the steady 
state in Fig. [7| as a function of the number of absorbed 
photons per chain. The frequency of the electric field is 
so varied again with interchain coupfings, ypi=0.3, 0.4, 
and 0.5 with C/p/V^i=0.9, that the energy of photons 
relative to that of the linear absorption peak is almost 
unchanged. The straight line shows the result previously 
obtained for isolated chains irrespective of the frequency, 
the amplitude, and the duration of the pulse.— ^ For finite 
interchain couplings, the final ionicity is again almost a 
linear function. With increasing interchain couplings, the 
final ionicity slightly deviates upward from the linear re- 
lation for large absorptions. It implies that they slightly 
enhance the cooperativity. Even for the largest couplings 
shown in the figure, however, the photoinduced ionic do- 
mains do not grow without energy supply. 



VI. SUMMARY 

Effects of interchain electron-electron interactions are 
studied on dynamics and nonlinear properties of phase 
transitions that are induced by intrachain charge-transfer 
photoexcitations, in a quasi-one-dimensional extended 
Hubbard model with alternating potentials and the 
electron-lattice coupling for mixed-stack charge-transfer 
complexes. Within the self-consistent mean-field approx- 
imation, we solve the time-dependent Schrodinger equa- 
tion for the electronic part and the classical equation of 
motion for the lattice part. 

In the ionic-to-neutral transition, the conversion frac- 
tion is a strongly nonlinear function of the photoexci- 
tation density for any interchain couplings. This co- 
operative property is maintained even if the electron- 
lattice interaction is turned offjiH. so that it is presum- 



ably due to electron-electron interactions that cause the 
ionic state basically a Mott insulator. All the electrons 
are so correlated in the ionic phase that any one of 
them cannot easily make a first move below the thresh- 
old. Because energies are accumulated by photoirra- 
diation of the ionic phase, each electron feels a force, 
but its motion is suppressed by electron-electron inter- 
actions. Once a metastable neutral domain is created 
above the threshold, it spontaneously grows leading to 
cooperative charge transfer. A similar collective charge- 
transport phenomenon is observed in field-effect transis- 
tors fabricated on organic single crystals of a quasi-one- 
dimensional Mott insulator j^ whose ambipolar charac- 
teristics are theoretically shown to be caused by balanc- 
ing the correlation effect in the bulk with the Schottky 
barrier effect at interfaces. ^^ 

For weak interchain electron-electron interactions, 
there is a range of the photoexcitation density that 
produces a mixture of neutral and ionic chains. This 
range is wide for small amplitude of the oscillating elec- 
tric field because the lattice-fluctuation-induced distri- 
bution of the threshold photoexcitation density is broad. 
For very large amplitude of the electric field, however, 
metastable domains are forced to grow simultaneously. 
Such coherence by the intense electric field survives only 
shortly. 

For strong interchain electron-electron interactions 
comparable to those in TTF-CA, the interchain corre- 
lation is strong during the transition. Neutral domains 
in nearby chains simultaneously grow even if their nucle- 
ation is delayed by reducing the amplitude of the electric 
field. Thus, the experimentally observed, coherent mo- 
tion of the macroscopic neutral-ionic domain boundary 



on a large time scaled is a consequence of the strong in- 
terchain couplings. 

In the ncutral-to-ionic transition, the conversion frac- 
tion is almost a linear function of the photoexcitation 
density. This property is maintained even when we pre- 
liminary added interchain elastic couplings to the present 
model: the energy scale of elastic interactions is so small 
compared with that of electron-electron interactions that 
elastic couplings would not largely modify the charge dy- 
namics at least on a short time scale. Such an uncooper- 
ative property is caused by uncorrelated electrons in the 
neutral phase, which move individually. The supplied 
energy is merely consumed to transfer electrons almost 
independently. The growth of metastable ionic domains 
is not spontaneous but always forced by the external field. 
Interchain electron-electron interactions make the func- 
tion slightly nonlinear, but the above situation is almost 
unchanged. This qualitative difference between the ionic- 
to-neutral and the neutral-to-ionic transitions is consis- 
tent with the experimental findings.^ 
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